Tracer Dispersion in a Self-Organized Critical System 
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We have studied experimentally transport properties in a 
slowly driven granular system which recently was shown to 
display self-organized criticality [Frette et al, Nature 379, 
49 (1996)]. Tracer particles were added to a pile and their 
transit times measured. The distribution of transit times is 
a constant with a crossover to a decaying power law. The 
average transport velocity decreases with system size. This is 
due to an increase in the active zone depth with system size. 
The relaxation processes generate coherently moving regions 
of grains mixed with convection. This picture is supported by 
considering transport in a ID cellular automaton modeling 
the experiment. 

PACS numbers: 05.40.+j, 64.60.Ht, 05.70.Jk, 05.70.Ln 



The avalanches that occur when grains are dropped 
onto a pile illustrate the spontaneous generation of com- 
plexity in simple dynamical systems Q. When grains 
are dropped onto a finite base, a pile builds up. How- 
ever, it cannot become infinitely high, and, eventually, 
the system settles in a stationary state where the out- 
flux over the edge of the base on average equals the in- 
flux. Intermittent flow of grains down the slope of the 
pile (small and large avalanches) maintain the system in 
this state. Bak, Tang, and Wiesenfeld constructed a 2D 
cellular automaton of a slowly driven dynamical system. 
They showed, that the "pile" spontaneously evolves, or 
self-organizes, into a state with avalanches of all sizes 
distributed according to a power law, that is, there is no 
internal system-specific scale. Because of the lack of any 
characteristic avalanche size, the system is referred to as 
critical H. 

It has been a longstanding question whether real 
granular systems display self-organized criticality (SOC) 
when slowly driven. Recently, however, an experiment on 
a quasi one-dimensional pile of rice has shown that the oc- 
currence of SOC depends on details in the grain-level dis- 
sipation mechanisms [|J . With nearly spherical grains, a 



characteristic avalanche size appeared, inconsistent with 
SOC. Only with sufficiently elongated grains, avalanches 
with a power-law distribution occurred. We focus in this 
Letter on the transport that results from avalanches in 
the system displaying SOC. The elongated rice grains 
could pack in a variety of ways, and each avalanche re- 
placed, locally or globally, one surface configuration with 
another. Thus a dynamically varying medium disorder 
(coupled to the relaxation processes) was generated. This 
is conceptually different from transport in media with a 
quenched disorder, see e.g. Refs. (3|^]. Furthermore, in 
SOC systems, a small perturbation may lead to arbitrar- 
ily large avalanches, and it is not clear at all, how this af- 
fects the transport properties. Thus it is quite surprising, 
that there are no experiments and only a few theoretical 
and numerical studies on transport in systems displaying 
SOCffJ. 

We have measured the transit times of colored tracer 
particles in the rice pile. Experimentally, we find that 
the distribution of transit times is essentially constant 
for small transit times T and decays as a power law for 
large T. A finite-size scaling analysis shows, that the av- 
erage transit time (T) oc L v where L is the system size 
and v — 1.5 ± 0.2. Thus the average velocity of tracer 
particles (V) oc L/(T) decreases with system size. This 
implies the existence of long-range correlations in the rice 
pile. Correlations in the sequence of transit times decay 
exponentially with a correlation time that increases with 
system size but have a crossover to a power-law decay. 
This indicates that the system has coherently moving re- 
gions, the size of which increases with system size. K\D 
cellular automaton model, related to the experiment, dis- 
plays SOC behavior and the distribution of transit times 
agrees qualitatively with the experimental findings. The 
scaling of the average transit time with system size is re- 
produced and hence the decrease of the average velocity 
with increasing system size. This phenomenon is related 
to an increase in the active zone depth as the system size 
is increased. Finally, the numerical and experimental re- 
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suits for the correlations in the transit times are very 
similar. 

The experimental system consisted of a rice pile con- 
fined between two 5 mm thick glass plates supported by 
15 mm thick 100 cm x 120 cm polymethylmethacrylate 
plates. Aluminum rods were inserted between the glass 
plates to form a vertical wall at one side and a vari- 
able base with length C of the quasi one-dimensional 
pile. The other side was open, allowing grains to fall 
off the pile. Grains of "Geisha Naturris" from Stabburet 
a.s. (Trollasen, Norway) with length 5 = 7.6 ± 0.9 mm 
and width 2.0 ± 0.1 mm were slowly fed into the gap be- 
tween the plates close to the vertical wall using an Ovrum 
Nibex 500 (Ovrum, Sweden) single seed machine. We 
used a plate separation of d = 6 mm and system sizes 
C = 15,30,60, and 85.7 cm, which, expressed in units of 
the grain length <5, are C/S = L = 20,39,79, and 113. 
The injection rate was 2-3 uncolored grains every 7.7 s 
or, on average, 20 grains/min. When the pile had reached 
the stationary state, color coded particles were added to 
the pile by hand and the transit times measured. A total 
number of 400 tracer particles for the two smaller systems 
and 800 for the two larger systems were used. Except for 
the smallest system, where a color coded tracer particle 
was injected every 2nd minute, the rate was 1 tracer par- 
ticle every 4th minute. The injection of uncolored grains 
was continued until all tracers had left the system. 

Figure 1 shows a part of the L = 113 system at a late 
stage of the experiment. Note that some of the tracer 
particles are close to the surface layer while others are 
buried quite deep in the pile. These tracers had to be 
released by large avalanches in order to move on. How- 
ever, large avalanches occurred with a very small prob- 
ability, and hence deeply buried particles tended to stay 
in the system for a very long time. Figure 2 is a full 
record of the experimental findings for the L = 113 rice 
pile. The projections of the horizontal lines onto the x- 
axis represent the time interval each tracer particle spent 
in the pile, that is, the transit time of the ith tracer 
T(i) = T out (i) - T in (i), where T in (i) and T out (i) denote 
the input and output time measured in units of addi- 
tions of uncolored grains (1 addition every 7.7 s), respec- 
tively. There is a huge variability in the transit times. 
The distribution functions of transit times P(T, L) for 
all system sizes are shown in Fig. 3(a). A data collapse 
for different system sizes L is obtained when plotting 
L 11 P(T,L) against the rescaled variable T/L" when us- 
ing v = 1.5 ± 0.2 and = 1.4 ± 0.2, see Fig. 3(b). Thus 
we can write 



P(T,L) = L- F(T/L V ), 



(1) 



where F is a scaling function and v a critical exponent 
expressing how the crossover transit time T c scales with 
system size. The scaling function F is of the form F(x) — 
constant for x < 1 and F(x) cx x~ a for a; > 1, where 



a = 2.4 ± 0.2. Since a > 1, it follows from the form of 
the scaling function F and the normalization constraint 
that (3 — v. The power-law tail does not contribute to the 
mean transit time (T) if a > 2, and (T) cx T c cx L v . In 
the experiments, the angle of repose was independent of 
system size. Thus the average velocity of tracer particles 
scales like (V) cx L/(T) oc L 1 ^ = L-°- 5±0 - 2 . It is quiet 
surprising, that the average velocity is not a constant but 
decreases with increasing system size. The (tracer) par- 
ticles have information on the system size! This can only 
happen if correlations exists throughout the system. In 
the statistically stationary state, the steady influx of par- 
ticles is balanced with a steady outflux. During the time 
At, each particle moves, on average, a distance At(V). 
Since the feeding rate was the same for all system sizes, 
the number of particles that crossed the outlet in this 
time interval was constant, that is, At (V) Xl = constant, 
where the active zone depth 
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h x {t) being the height at position x at time t and ( • } t 
denotes the temporal average. This implies that the av- 
erage velocity is inversely proportional to the active zone 
depth, (V) cx 1/Xl- An increase in Al with system size 
would be consistent with the experimental finding of a 
decreasing average velocity with increasing system size. 

Further insight into the correlations that exist in this 
dynamical state can be extracted from the data. We 
notice that, once in a while, many tracers dropped out 
of the system at the same time since they were part of a 
large avalanche that reached the rim of the pile. The cor- 
responding "steps" in the diagrams, like the one shown 
in Fig. 2, have many different sizes and are interwoven in 
a complex way. The correlation function 



C (t) = (I(i,i + T)), 



where 
J(M + t) 



1 if \T out (i) 
otherwise 



T out (i + T)\ < St 



(3) 



(4) 



is the indicator function of simultaneously drop out and 
( • } denotes average over all i, gives the probability that 
tracer i and i + t left the system simultaneously. Us- 
ing St = 10, this quantity is found to behave as c(r) cx 
exp(— t/t c (L)) for small r, where r is the time difference 
in units of 4 minutes and t c (L) = 5.8± 1.3, 8.7 ±2, 19 ±4, 
and 34 ± 4 is the correlation time for system sizes L = 
20,39,79, and 113, respectively. In Fig. 2, we see that 
t c (L = 113) = 34 corresponds to the average size of the 
steps. Since the injection rate was 20 grains/min, the 
characteristic number of grains spanned by these corre- 
lated sequences of tracers were 476, 705, 1539, and 2754. 
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The solid-like motion of domains described in Ref. ||] 
is probably one aspect of this coherent motion. Using 
Xl oc volume/L and disregarding the smallest system, 
we obtain Al oc £°- 25 ±°- 2 5 consistent within error bars 
with the result above. However, there is significant dis- 
persion. Grains are being transferred between different 
coherently moving areas, see Fig. 2. 

Inspired by the experiments, we considered a refined 
version of a simple ID cellular automaton studied in Ref. 
p|. In a system of size L, an integer variable h x gives the 
height of the pile at site x. The local slope z x at site x is 
given by z x = h x — h x+ i, and we impose /il+i = 0. The 
addition of a grain at the wall increases the slope by one 
at x = 1, that is, z\ — ► Z\ + 1. We proceed by dropping 
grains at the wall until the slope z\ exceeds a critical 
value, z\ > zl, then the site topples by transferring one 
grain to its neighboring site on the right, x = 2. If z x > 
z^, this site topples in turn according to 

!*"iV 2 +1 (5) 

z x ±i — ► z x ±i + 1 

(unless at the rightmost site where the grains fall off the 
pile) generating an avalanche. During the avalanche, no 
grains are added to the pile. Thus the two time scales in- 
volved in the dynamic evolution of the pile are separated. 
The injection rate of grains is low compared to the du- 
ration of the relaxation processes. The avalanche stops 
when the system reaches a stable state with z x < z^ \fx 
and grains are added at the wall until a new avalanche is 
initiated and so on. The critical slopes z x are dynamical 
variables chosen randomly to be 1 or 2 every time site x 
has toppled. This is a simple way to model the changes 
in the local slopes observed in the rice pile experiment. 
Thus the model differs from the trivial II? Bak, Tang, 
Wiesenfeld model where z c x = 1 is a constant and where 
grains are added on randomly chosen sites [[l]. The ran- 
domness in the BTW model is external. In our model, 
the randomness is internal and inherent in the dynamics. 
Starting with, say, z x = and z% = 1 Vx, the system 
reaches a stationary state where the avalanche sizes are 
power-law distributed with an exponent of —1.55 ± 0.10 
and a cutoff in the power-law distribution that scales with 
system size as £ 2 - 25±al °. 

We have measured the transit times (in units of added 
grains) of the added particles after the pile has reached 
the stationary critical state. Using Eq. (Q), a reasonable 
data collapse is obtained with v = 1.30 ± 0.10 and (3 = 
1.35±0.10 as displayed in Fig. 4. We find a = 2.22±0.10, 
that is, the results agree well with the experimental find- 
ings. The average depth of tracer particles during the 
transport through the system as a function of the aver- 
age transit time is displayed as an inset of Fig. 4, and it 
shows, that the power-law tail in the distribution of tran- 
sit times arises from tracer particles that become deeply 
buried in the pile, while tracer particles that stayed close 



to the surface of the pile during the transport through 
the system mainly contributed to the constant part. The 
average velocity (V) oc L/(T) oc L 1 ^ = L~°- 30±0 W , 
in fairly good agreement with the experimental result. 
Furthermore, in the model, we are able to measure di- 
rectly the scaling of the active zone depth with system 
size. We find A^ oc £°- 25±0 - 10 in agreement with the re- 
lation (V) oc 1/Al- Alternatively, a third measure for 
the active zone depth is available from the inset in Fig. 
4, giving Al oc depth oc £°- 3±010 . Moreover, correla- 
tion functions calculated from Eqs. (3) and (4) give an 
exponential decay for short times with a crossover to a 
decaying power law for large times. Adding grains at the 
wall and allowing for dynamical critical slopes, we have 
seen similar results for the nonlocal limited model in Ref. 
. Thus the behavior seems to be universal. 
In conclusion, this new direction of research sheds light 
upon the dynamics of SOC systems in general and gran- 
ular systems in particular. We find that transport in a 
self-organized critical granular medium is characterized 
by an average grain velocity that approaches zero when 
the system size increases. The dynamics of the system 
displays comprehensive correlations. These experimen- 
tal findings agree well with the behavior seen in simple 
one-dimensional computer models of the self-organized 
critical pile. 
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FIG. 1. A close up photograph of the rice pile with size 
L = 113. The plate separation/grain length ratio d/S = 0.79 
and most of the grains are aligned along the flow direction. 
Each tracer particle is uniquely color coded; thus it is possible 
to measure individual transit times. This picture was taken 
after 638 out of 800 tracer particles had been added. The 
region shown contains 5 of a total of 73 tracer particles that 
were inside the pile at this time. 

FIG. 2. A record of the tracer experiment in a pile of size 
L = 113 where a total number of 800 tracer particles were 
added, one every 4th minute. The tilted line connects the 
injection times for all the tracers. The transit time for each 
tracer particle is represented by the length of a horizontal 
line whose projection onto the x-axis of the left (right) end- 
point equals the time the particle entered T; n (left T ou t) the 
system. The transit time is measured in units of the num- 
ber of injections of uncolored grains (no. additions), 1 addi- 
tion every 7.7 s. Note the large variation in the transit times 
T — T ou t — Ti n and that, repeatedly, many tracers left the 
system at the same time. The horizontal arrow indicates the 
time at which the photograph in Fig.l was taken. Further- 
more, the correlation time t c (see text related to Eq. (^)) is 
indicated. The inset is the part of the full record marked by 
dashed lines and shows that the different sizes of steps are 
interwoven. 

FIG. 3. (a) The experimental results for the normalized 
distribution of transit times in piles with sizes L = 20, 39, 79, 
and 113. The data have been averaged over exponentially in- 
creasing bins with base 2 in order to reduce the fluctuations in 
the statistics due to the relatively small number of tracer par- 
ticles. The functions are essentially constant for small T and 
have a decaying power-law tail with a slope of a = 2.4 ± 0.2. 
These large transit times correspond to tracer particles which, 
during the transport through the systems, become deeply em- 
bedded in the pile, (b) Disregarding the smallest system and 
using Eq. ([!]), a reasonable data collapse of the three largest 
systems is obtained with v = 1.5 ± 0.2 and j3 = 1.4 ± 0.2. 

FIG. 4. A finite-size scaling plot using Eq. (^) with 
v = 1.30±0.10 and (3 = 1.35±0.10 of the normalized distribu- 
tion of transit times in the numerical model with system sizes 
L — 25, 100, 400, and 1600. The statistics shown correspond 
to 10 7 tracer particles (10 6 for L = 1600), and the data have 
been averaged over exponentially increasing bins with base 
1.1. The functions are constant for small transit times and 
decay as power laws with a slope of a = 2.22 ± 0.10. The 
inset shows that a data collapse of the correlation between 
the transit times and depth of the tracer particles can be ob- 
tained using a finite-size scaling plot analogous to Eq. (|l|), 
with v' = 1.20 ± 0.10 and fi' = -0.30 ± 0.10. 



